Topological transport in Dirac electronic systems: A concise review
Song Hua-Ding1, Sheng Dian1, Wang An-Qi2, Li Jin-Guang1, Yu Da-Peng1, 3, 4, Liao Zhi-Min1, 3, †
State Key Laboratory for Mesoscopic Physics, School of Physics, Peking University, Beijing 100871, China
Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
Collaborative Innovation Center of Quantum Matter, Beijing 100871, China
Department of Physics, South University of Science and Technology of China, Shenzhen 518055, China

 

† Corresponding author. E-mail: liaozm@pku.edu.cn

Abstract

Various novel physical properties have emerged in Dirac electronic systems, especially the topological characters protected by symmetry. Current studies on these systems have been greatly promoted by the intuitive concepts of Berry phase and Berry curvature, which provide precise definitions of the topological phases. In this topical review, transport properties of topological insulator ( , topological Dirac semimetal ( ), and topological insulator-graphene heterojunction are presented and discussed. Perspectives about transport properties of two-dimensional topological nontrivial systems, including topological edge transport, topological valley transport, and topological Weyl semimetals, are provided.

1. Introduction

Inspired by the notion of “topological order”, the exploration for novel classifications of matter phases becomes an intriguing theme in condensed-matter physics. Over three decades, the discovery of the quantum Hall effect (QHE) enlightened investigations on topological phases of matters.[16] Assisted by the powerful concepts of Berry phase and Berry curvature, precise definitions of topological phases are achieved. Berry phase is the phase generation during an adiabatic evolution process with steady external parameters producing a loop in the parameter space.[1,2] The Berry phase can be written as the integral of the Berry curvature.[2] Moreover, the integrals of the Berry curvature over closed surfaces are always employed to identify topological nontrivial phases, known as quantized Chern numbers.[2] Hence, theoretical predictions for nontrivial matters emerged endlessly for more than a decade, and novel quantum systems were continuously revealed. Such topological nontrivial systems exhibit fascinating properties, especially in transport experiments.

The QHE was discovered by Klaus von Klitzing et al., where the Hall conductivity of a two-dimensional (2D) system was found to be exactly quantized under a strong magnetic field.[79] Thouless et al. firstly employed topological considerations to explain the QHE.[5] From their milestone work, researchers began to realize that the electronic topological structure of quantum materials may be directly related to the transport behaviors. Various materials with novel topological phase have already been revealed. Graphene, the first 2D material containing Dirac fermions, has been proved to possess a Berry phase of π, leading to half-integer QHE observed in transport measurements.[1016] Topological insulators (TIs) like and , where strong spin-orbit coupling creates topological insulating electronic phases and non-zero integer Chern numbers, are very potential platforms for spintronics and next-generation information devices.[1733] Gapless surface states protected by time-reversal symmetry exist on the surface of these materials, and less-dissipative transport can be achieved by the surface states, which will lead to an information highway in the future. Similar to the graphene system, the surface states of TIs also behave as 2D Dirac fermions with π Berry phase.[1821] Thus, Hall conductivity of the surface states is half-integer quantized as well. However, the parallel connection of the top and bottom surfaces may necessarily carry an integer QHE in real measurements.[23] Rui Yu et al. proposed that quantum anomalous Hall effect (QAHE) can be realized in a magnetized TI system,[25] which was confirmed by Cui-Zu Chang et al.[26,27] Recently, topological semimetals became another focus of condensed-matter physics. The topological Dirac semimetals (for example, Bi and ) are natural three-dimensional counterpart of graphene, with bulk linear-dispersion Dirac fermions and Fermi arcs surface state.[3439] The unique chiral character of the Dirac fermions in topological Dirac semimetals generates abundant physical entities, such as chiral anomaly induced charge pumping behavior.[39] Novel transport properties, like anomalous negative magnetoresistance and surface state assisted Aharonov–Bohm oscillations,[40,41] were successively discovered.

In this topical review, we generally introduce the experimental investigations about transport behaviors in topological nontrivial systems, mainly based on our own works. We have reported abundant phenomena originating from surface states in topological insulators, topological Dirac semimetal, and monolayer graphene. In Section 2, we demonstrate the enhanced photothermoelectric effect by TI helical surface states[42] and the quantum oscillations from TI nanoribbon sidewalls.[43] In Section 3, our investigations on the topological Dirac semimetal ), including negative magnetoresistance induced by chiral anomaly,[40] Aharonov–Bohm oscillations in nanowires, and two-carrier transport in nanoplates are summarized.[41,44] In Section 4, we exhibit the novel tunneling behavior in the vertical heterojunction formed by graphene and TI.[45] In Section 5, we will review other significant theoretical and experimental progresses on the topological edge state transport, the valley transport, and the topological Weyl semimetals. Finally, we give perspectives and conclusions in Section 6.

2. Transport properties of topological insulator
2.1. Enhanced photothermoelectric effect of 2D topological surface states

Here we review the photothermoelectric effect in topological insulator nanoribbons. , as a strong three-dimensional topological insulator, has a bulk energy gap of ∼ 0.3 eV and helically conducting chiral surface states.[22,24] The 2D surface states could be excited to be spin-polarized by means of circularly polarized light. Due to spin-momentum locking in topological insulators, the spin-polarized surface states have oriented motions, which could be accelerated by the temperature gradient, resulting in an enhanced photothermoelectric effect.[46,47]

The nanoribbons were synthesized via the chemical vapor deposition (CVD) method. To better characterize the crystalline quality and study the layered microstructure of the nanoribbons, the as-synthesized nanoribbons were prepared to be cross-sectional on a substrate (Fig. 1(a)). Measured from the cross section along the [100] zone axis, the diffraction pattern indicates that the nanoribbon is single crystalline (Fig. 1(b)) with the basic unit of quintuple layer of –Se–Bi–Se–Bi–Se– sequentially (Fig. 1(c)). The ARPES spectrum acquired at room temperature confirms the existence of TI surface states (Fig. 1(d)).

Fig. 1. (color online) (a) Cross-sectional HAADF-STEM image of a nanoribbon. (b) Diffraction pattern acquired from the cross-section along the [100] zone axis. (c) High resolution HAADF-STEM image of the nanoribbon surface, where the quintuple layer layer is indicated. (Bi represented by red dots; Se represented by green dots.) (d) Room-temperature ARPES spectrum for a nanoribbon.[42]

Individual nanoribbons were fabricated to contact with Cr/Au electrodes on a Si substrate with 300 nm layer via electron beam lithography techniques, as shown in Fig. 2(a). To study the surface state enhanced photothermoelectric effect in the nanoribbons, the voltage response to circularly polarized light illumination was measured. Right-circular polarization (RCP), linear polarization (LP), and left-circular polarization (LCP) of the incident laser were achieved via tuning the λ/4 waveplate. The incident laser was along the yz plane, where θ, the relative orientation of the incident light and the z axis, was chosen to be 0° and 30° (Fig. 2(b)). Considering that the most notable voltage response was detected when illuminating near the electrodes, the incident laser in our experiment was illuminated near the measurement electrodes. At θ = 0°, the voltage vs. time curves of RCP, LP, LCP show a little difference, as a result of that the vertically incident laser cannot induce in-plane surface states to be spin polarized, thus no additional surface state related photocurrent can be excited.

Fig. 2. (color online) (a) Optical image of a typical nanoribbon PTE device. (b) Sketch of the measurement set up. (c) The photo-voltage response measured when the laser illumination on was near the positive electrode. (d) The photo-voltage response measured when the laser illumination on was near the ground electrode. (e) The physical mechanism for the light polarization-selective transition.[42]

At θ = 30°, it is found that the RCP and LCP light have different effects on the photothermoelectric output. When the laser spot is near the electrode connected to the voltmeter positive terminal (the positive electrode), the RCP light induces an obviously enhanced voltage signal, the voltages induced by LCP and LP are similar (Fig. 2(c)). However, when the laser is illuminated near the voltmeter ground terminal (the negative electrode), the LCP induced photothermoelectric effect is notable, while the voltage response of RCP is almost the same as that of LP (Fig. 2(d)).

The phenomena above are dominated by the light polarization selective rule according to the helical character of the carriers (Fig. 2(e)). We define the electrons in the Dirac cone with a positive slope as spin-up, while those in the branch with a negative slope correspond to spin-down states. The quantum numbers of the orbital angular momentum of each branch are +1 and −1, respectively.[4851] A photon can selectively excite the surface states to the conduction band states (l = 0) by transferring angular momentum.[48,51] If spin-polarized surface electrons have the same direction as the electrons driven by the temperature gradient, the generated voltage would be enhanced. Conversely, if the additional surface electrons have the opposite motion direction with the electrons driven by the temperature gradient, there would be little impact on the total generated voltage.

When the laser illumination on is near the positive electrode, the σ– photons (RCP) only excite the surface electrons with to conduction bands (l = 0), while the hopping of electrons is forbidden, resulting in additional spin-up surface states. Considering the spin-momentum locking in , the spin-up surface electrons have a translational direction along the +x axis,[52] which is parallel to the direction of the temperature gradient. Naturally, additional surface electrons would be accelerated and enhance the photothermoelectric effect. Contrarily, the LCP would hardly enhance photovoltage as a result of the anti-parallel relationship between the momentum of spin-down surface electrons and the temperature gradient. A similar analysis can be made for the laser illumination position near the ground electrode, in this case, the LCP would induce significant enhancement of the photothermoelectric effect.

2.2. Side-wall surface states associated quantum oscillations

Three-dimensional topological insulators have been demonstrated to possess rich physical connotations.[1724] is always considered as a perfect platform to explore the novel properties of 2D Dirac fermions on its surface. However, even containing a large bulk energy gap, the conductance contributed from the bulk is still notable due to Se vacancies in making identification of the surface states by transport measurements a great challenge. In this section we review our investigations on prominent SdH oscillations originating from the surface Dirac fermions in the sidewalls of nanoplates. Importantly, the SdH oscillations from the sidewalls appear with a dramatically weakened magnetoresistance background, offering a direct path to detect the 2D surface states when the coexistence of bulk and surface conductions is inevitable.

The nanoplates were prepared by the vapor-liquid-solid mechanism, using gold as the catalyst. Electrodes of Cr/Au (10 nm/130 nm) were deposited by electron beam evaporation for transport measurements (Fig. 3(b)). When applying an in-plane magnetic field, obvious SdH oscillations can still be observed on a smoothed magnetoresistance background (Fig. 3(c)), the oscillation features are gradually smeared out by increasing temperature (over 20 K). Such in-plane field induced SdH oscillations have rarely been reported before, and apparently, this phenomenon does not originate from the top/bottom surfaces. These results indicate that the most possible origination is the surface states on two sidewalls.

Fig. 3. (color online) (a) Schematic of the in-plane magnetic field referring to the sidewall (Se: green spheres; Bi: purple spheres).(b) Optical image of a typical nanoplate device and a schematic of the orientation of the applied magnetic field. The length between the two inner longitudinal voltage leads, and the width of the two Hall voltage leads of the nanoplate device are 1.48 μm and 2.26 μm, respectively. (c) Plot of longitudinal resistance versus under in-plane magnetic field at different temperatures. (d) Landau-level fan diagram for the Shubnikov–de Haas oscillation in at the lowest temperature 1.3 K.[43]

Considering the Onsager semiclassical quantization relation, a linear relationship can be acquired: , where indexes the N-th minimum in conductance, K is a constant, and is related to the Berry phase by . Based on such a relation, the Landau level fan diagram can be established (Fig. 3(d)). The extrapolated Landau index at extreme field limit ( ) is about 0.5, acquired from the Landau level fan diagram, indicating an additional Berry phase π in the system. This is direct evidence of novel 2D Dirac fermions. We substantially demonstrated that the oscillation behavior derived from Landau level formation in 2D Dirac fermions. Since the applied in-plane field cannot induce such formation in the top/bottom surface, the sidewalls are the only possible origin for the quantum oscillation we observed.

Angle dependent magneto-transport measurements show that the positive magnetoresistance background, which is the contribution from three-dimensional bulk transport, remarkably decreases as the out-of-plane field component gradually disappears (Fig. 4(a)). To avoid the complexity that top/bottom and sidewalls simultaneously take part in transport, we measured the transversal resistance (2 K), which should not be affected by the sidewall conductance (Fig. 4(b)). Clear oscillation features have been observed (Fig. 4(c)). We derived the Landau level fan diagrams from the results in Fig. 4(b). We found that the derived slopes can be well fitted by 1/cos θ (Fig. 4(d)). Therefore we can confirm that the observed SdH oscillations avoid influences from three-dimensional bulk transport, and the oscillations under in-plane field should originate from the sidewalls. Moreover, such sidewall associated SdH oscillations successfully avoid disturbances from bulk transport, providing an easy method to detect surface states in topological insulators.

Fig. 4. (color online) (a) Longitudinal resistance and (b) Hall resistance as a function of the applied magnetic field B under different orientations described by the angle θ at a temperature of 2 K. (c) Oscillations of the Hall resistance at 2 K and at different angle θ. (d) Fitted slope of the Landau level fan diagram derived from as a function of the angle θ; the solid line gives a good fitting with 1/cos θ.[43]
3. Transport properties of topological Dirac semimetal
3.1. Giant negative magnetoresistance induced by the chiral anomaly

is one kind of the so-called three-dimensional Dirac semimetals, which are also described as three-dimensional analogues of graphene. These materials are unusual quantum materials containing massless Dirac fermions.[19] By breaking the time-reversal symmetry or spatial inversion symmetry, the Dirac semimetal is believed to transform into a Weyl semimetal with chiral anomaly.[3439,53] In single crystal nanoplates, we find large negative magnetoresistance induced by chiral anomaly. Moreover, the negative magnetoresistance is closely linked with the carrier density, which is tunable by gate voltage and temperature. Our finding is helpful to understand Weyl fermions in Dirac semimetals.

The nanoplates were prepared by the chemical vapor deposition method. The exposed surface of the nanoplate samples was the (112) plane, identified by TEM results (Fig. 5(a)). The nanoplates possess a thickness of several hundred nanometers (Fig. 5(b)). The samples were then fabricated to contact with Au electrodes on an oxide Si substrate, serving as the back gate. Four-probe measurements were accepted in the experiments (Fig. 5(c) inset). The temperature dependence of resistivity of the samples demonstrates semiconducting-like behavior, which is different from the metallic behavior in bulk crystal (Fig. 5(c)). This is mainly attributed to the deviation of carrier density. The sample has a low carrier density and its Fermi level is close to Dirac points, where the holes in valence bands can be easily thermally activated to conduction bands in high temperature. So the resistivity of increases as the temperature decreases due to the reduced thermal activation. After decreasing to a critical low temperature, the thermal energy is not enough for holes to activate across the gap, leading to the metallic behavior. Unlike previous nanowire, the bulks usually maintain high carrier density and the Fermi level is always above Dirac points, where the metallic pT mechanism is dominant. As shown in Fig. 5(d), nanoplates exhibit a notable negative MR with in-plane B, even at room temperature. The negative MR also emerges in nanowire samples with . The negative MR is mainly induced by the chiral anomaly effect. , as a three-dimensional Dirac semimetal, its Dirac point is composed of two overlapping Weyl nodes with opposite chiralities (left handed and right handed) (Fig. 6(a)), which would be separated in the magnetic field.[54,55] Applying an extra parallel electric field, the two kinds of Weyl fermions would have unequal chemical potentials ( ) (Fig. 6(b)). In such a situation, the continuity equation of the Weyl nodes is

(1)
In addition to the typical band dispersion contribution, an anomalous term due to the Berry curvature emerges, the extra term is always transverse to the electric field, leading to a Hall current.[2] For the system with inversion symmetry breaking, the Berry curvatures in different valleys have opposite signs under the protection of time-reversal symmetry, thus non-zero valley conductivity can exist in the transverse direction, a net valley degree accumulation is expected on both transverse sides. This phenomenon is called the “valley Hall effect”. To experimentally realize the VHE, suppressed valley relaxation must be achieved. The spontaneous inversion symmetry breaking and strong spin–orbit coupling make an ideal platform to realize VHE.[90] Mak et al. firstly reported the experimental observation of VHE in monolayer .[91] Circularly polarized light was used to create a chemical potential imbalance in the two valleys of , inverting the valley signal to an observable electric signal; an anomalous Hall signal induced by VHE was substantially detected.

5.2.2. The valley transport in graphene

Theoretical analysis has shown that the graphene with inversion symmetry breaking can also exhibit the VHE.[92] However, pristine graphene (single layer or bilayer) is under inversion symmetry protection. The experimental investigations on the valley transport in a graphene system have been achieved recently, with the help of improvement on device fabrication.

In 2014, Gorbachev et al. reported a large non-local signal detected in single layer graphene with underlying BN substrate.[93] The inversion symmetry of graphene is believed to be broken as the consequence of A/B sublattice global asymmetry induced by the substrate potential. The inversion symmetry breaking leads to a non-zero Berry curvature distribution in the Brillouin zone. When the Fermi level is tuned to the energy level where non-zero Berry curvature exists, a large non-local signal can be detected from the transport measurements. The anomalously obvious non-local signal is due to the valley Hall effect and the inverse VHE, the latter is the phenomenon where pure valley current gives rise to charge current.

The novel work of Gorbachev et al. not only firstly reported the experimental realization of graphene valley transport, but also provided an efficient method to detect the VHE. In 2015, the valley transport in bilayer graphene was discovered by Sui et al. and Shimazaki et al. The two groups both utilized a perpendicular electric field tuned by a gate voltage to break the inversion symmetry of graphene, a strong non-local signal was detected, which was attributed to the VHE and the inverse VHE.[94,95] Their investigations show that electrically driven inversion symmetry breaking is achievable, shedding light on future valleytronics. Furthermore, the AB/BA stacking domain wall in bilayer graphene has also been reported to support a non-diffusional valley current.[96] The bilayer graphene system with broken inversion symmetry may possess a non-zero “valley” Chern number, leading to a topological valley edge state existing along the domain wall.

5.3. Topological Weyl semimetals

Based on the quantum field theory, the fermions of our universe are divided into three kinds: the Dirac fermions, the Weyl fermions, and the Majorana fermions.[97] Interestingly, the condensed-matter systems host quasiparticles possessing the characters of these high-energy physics concepts. The most typical example of Dirac semimetals is graphene, which has a Dirac point protected by time-reversal symmetry and inversion symmetry.[1016] Extending the 2D linear energy–momentum dispersion relation to 3D momentum space, the 3D Dirac semimetals were recently discovered. Moreover, the Dirac point in 3D Dirac semimetals is composed of two Weyl nodes. The topological Weyl semimetals are the systems where low-energy excitations obey the Weyl equations.[57,98100] Compared to the Dirac points (four-fold degeneracy), the Weyl nodes are only two-fold degenerate. Thus the Weyl node can only exist in systems with spontaneous symmetry-breaking, because in the systems possessing both the time-reversal and inversion symmetries, the band crossings are always four-fold degenerate.

The topology of the Weyl nodes is also defined by the Berry curvature,[98100] an analogy of the magnetic flux in momentum space. Generally speaking, the Weyl nodes can be considered as the “magnetic monopoles” in momentum space, which are the origins or ends of the Berry curvature. The chirality (+1 or −1) of a Weyl node is similar to the “magnetic charge”, reflecting the Berry curvature distribution embracing the Weyl node. The Weyl nodes always come in pairs because the Chern number of the Weyl semimetal must be zero, otherwise it becomes a topological insulator.[2]

There are a large number of unusual physical features associated with Weyl semimetal, including open Fermi arcs in the surface and various exotic transport phenomena induced by the chiral anomaly. The open Fermi arcs connect the projections of the bulk Weyl nodes with opposite chiralities onto the surface, and the Fermi arcs can be understood as the surface states of Weyl semimetals.[98] Under the presence of external electric and magnetic fields, the chiral charge of one Weyl node is not conserved, which is the so-called chiral anomaly. The charge pumping effect between Weyl nodes with opposite chiralities gives rise to large negative magnetoresistance.[40,5457,97,98]

It is worth noting that the Weyl semimetals discussed above have zero density of state at the zero energy, and they obey the Lorentz invariance in quantum field theory. This kind of Weyl semimetals have been predicted to exist in several symmetry-breaking materials, including TaAs, NbP, and other systems.[99,100] By the ARPES technique, the existence of Weyl fermions in TaAs, NbAs has been experimentally confirmed,[101105] and chiral-anomaly-induced negative magnetoresistance was also detected in TaAs recently.[106]

Recently, a new type of Weyl semimetal has been predicted in the system.[107] In this type of Weyl semimetal, non-zero density of state at zero energy exists, and such system seems to break through the limit of Lorentz invariance, which is referred as a type-2 Weyl semimetal. The remarkable distinction between the Fermi surfaces of the two types of Weyl semimetals leads to notable differences in the thermodynamics and magneto-transport. The type-2 Weyl semimetal exhibits the chiral anomaly effect with strong anisotropy.

6. Conclusion and perspectives

Topological nontrivial systems provide splendid platforms for uncovering novel physics in mesoscopic scale, and also for practical applications. As a powerful method of physical characterization, transport measurements make significant contributions in exploring these miraculous systems. This topical review summarizes our experimental investigations on transport behaviors in topological insulator ( , topological Dirac semimetal ( ), and TI-graphene heterojunction. We demonstrate fantastic quantum effects of these systems, exhibiting their great potential in realizing next-generation electronic and spintronic devices. We also review significant breakthroughs from other groups, to depict a general status of current research on topological transport. Although great triumphs have been achieved, future investigations are still in eagerly awaited.

The past ten years have witnessed great achievements on the graphene system. Moreover, this intriguing 2D material still possesses unpredictable academic and application prospects. Researchers begin to enlarge the functionalization of graphene using graphene-based heterostructures by the proximity effect. Other than the inversion symmetry breaking and enhancement of spin–orbit coupling discussed above, a large exchange field has been introduced into graphene by a magnetic insulator EuS.[108] We believe experimental realizations of QSHE and QAHE in graphene-based heterostructures will come true in the near future. The physical contents and application possibilities are obviously diversified by investigations on the graphene-based heterostructures.

The proximity effect also shows extraordinary talents in other two-dimensional systems. The superconducting proximity effect is believed to induce Majorana fermions in a topological nontrivial system,[109] and transform the hybrid system into a “topological superconductor”. The controllable realization of such fantastic fermions will start the revolution of quantum computing. These features lead the superconducting proximity effect to become a research focus in recent years. Researchers have successfully detected the signature of Majorana fermions in topological superconductor nanowires,[110,111] and we notice that a recent report has provided evidence of the Majorana zero mode in the topological superconductor.[112] The superconductivity in Dirac semimetal was also achieved by a point contact.[113] The topological superconductor systems will remain intriguing for a long time, and there is still an inestimable amount of physical content for us to discover.

Many other novel 2D materials are waiting to be revealed. The black phosphorus, which has both band gap and high mobility, was discovered in recent years.[114117] Researchers have inquired into its application prospect,[114,115] and the QHE has been realized in black phosphorus.[116,117] Other new systems, like type-2 Weyl semimetal, are explored continuously. The variety of 2D materials affords us incessant inspirations, and the combination of 2D materials and other systems (for example, 2D material with superconductors or ferromagnets) will actually maximize such variety.

On the other hand, the improvement of device fabrication techniques provides substantial foundation for intensive works. The MBE method has been developed to grow high-quality 2D materials and their Van der Waals heterostructures. The interference of impurity potential can be greatly suppressed by the device encapsulation technology (for example, BN encapsulation),[118] leading to an ultraclean device environment. By the help of technological advances, more and more physical phenomena have been discovered experimentally.

We hope the readers of Chinese Physics B may acquire inspirations from our review, and our investigations over these years may shed light on future explorations.

Reference
[1] Berry M V 1984 Proc. R. Soc. London Ser. A 392 45
[2] Xiao D Chang M C Niu Q 2010 Rev. Mod. Phys. 82 1959
[3] Kohmoto M 1985 Ann. Phys. 160 343
[4] Simon B 1983 Phys. Rev. Lett. 51 2167
[5] Thouless D Kohmoto M Nightingale M Den Nijs M 1982 Phys. Rev. Lett. 49 405
[6] Wilczek F Shapere A 1989 Geometric Phases in Physics Singapore World Scientific
[7] Klitzing K V Dorda G Pepper M 1980 Phys. Rev. Lett. 45 494
[8] Mendez E E Esaki L Chang L L 1985 Phys. Rev. Lett. 55 2216
[9] Tsui D C Stormer H L Gossard A C 1982 Phys. Rev. Lett. 48 1559
[10] Castro Neto A H Peres N M R Novoselov K S Geim A K 2009 Rev. Mod. Phys. 81 109
[11] Das Sarma S Adam S Hwang E H Rossi E 2011 Rev. Mod. Phys. 83 407
[12] Geim A K Novoselov K S 2007 Nat. Mater. 6 183
[13] Novoselov K S Geim A K Morozov S Jiang D Zhang Y Dubonos S Grigorieva I Firsov A 2004 Science 306 666
[14] Novoselov K S Geim A K Morozov S V Jiang D Katsnelson M I Grigorieva I V Dubonos S V Firsov A A 2005 Nature 438 197
[15] Novoselov K S Jiang D Schedin F Booth T J Khotkevich V V Morozov S V Geim A K 2005 Proc. Natl. Acad. Sci. U S A 102 10451
[16] Zhang Y Tan Y W Stormer H L Kim P 2005 Nature 438 201
[17] Chen Y Analytis J Chu J H Liu Z Mo S K Qi X L Zhang H Lu D Dai X Fang Z 2009 Science 325 178
[18] Fu L Kane C L 2007 Phys. Rev. B 76 045302
[19] Hasan M Z Kane C L 2010 Rev. Mod. Phys. 82 3045
[20] Wang J DaSilva A M Chang C Z He K Jain J Samarth N Ma X C Xue Q K Chan M H 2011 Phys. Rev. B 83 245438
[21] Qi X L Zhang S C 2011 Rev. Mod. Phys. 83 1057
[22] Xia Y Qian D Hsieh D Wray L Pal A Lin H Bansil A Grauer D Hor Y Cava R 2009 Nat. Phy. 5 398
[23] Xu Y Miotkowski I Liu C Tian J Nam H Alidoust N Hu J Shih C K Hasan M Z Chen Y P 2014 Nat. Phy. 10 956
[24] Zhang H Liu C X Qi X L Dai X Fang Z Zhang S C 2009 Nat. Phy. 5 438
[25] Yu R Zhang W Zhang H J Zhang S C Dai X Fang Z 2010 Science 329 61
[26] Chang C Z Zhang J Feng X Shen J Zhang Z Guo M Li K Ou Y Wei P Wang L L Ji Z Q Feng Y Ji S Chen X Jia J Dai X Fang Z Zhang S C He K Wang Y Lu L Ma X C Xue Q K 2013 Science 340 167
[27] Chang C Z Zhao W Kim D Y Zhang H Assaf B A Heiman D Zhang S C Liu C Chan M H Moodera J S 2015 Nat. Mater. 14 473
[28] Xie M H Guo X Xu Z J Ho W K 2013 Chin. Phys. B 22 068101
[29] Gao C L Qian D Liu C H Jia J F Liu F 2013 Chin. Phys. B 22 067304
[30] Wang J Zhu B F 2013 Chin. Phys. B 22 067301
[31] Liu Y Ma Z Zhao Y F Meenakshi S Wang J 2013 Chin. Phys. B 22 067302
[32] Sheng L Li H C Yang Y Y Sheng D N Xing D Y 2013 Chin. Phys. B 22 067201
[33] He K Ma X C Chen X Lu L Wang Y Y Xue Q K 2013 Chin. Phys. B 22 067305
[34] Liu Z Jiang J Zhou B Wang Z Zhang Y Weng H Prabhakaran D Mo S Peng H Dudin P 2014 Nat. Mater. 13 677
[35] Liu Z Zhou B Zhang Y Wang Z Weng H Prabhakaran D Mo S K Shen Z Fang Z Dai X 2014 Science 343 864
[36] Neupane M Xu S Y Sankar R Alidoust N Bian G Liu C Belopolski I Chang T R Jeng H T Lin H Bansil A Chou F Hasan M Z 2014 Nat. Commun. 5 3786
[37] Wang Z Sun Y Chen X Q Franchini C Xu G Weng H Dai X Fang Z 2012 Phys. Rev. B 85 195320
[38] Xiong J Kushwaha S K Liang T Krizan J W Hirschberger M Wang W Cava R J Ong N P 2015 Science 350 413
[39] Parameswaran S Grover T Abanin D Pesin D Vishwanath A 2014 Phy. Rev. X 4 031035
[40] Li C Z Wang L X Liu H Wang J Liao Z M Yu D P 2015 Nat. Commun. 6 10137
[41] Wang L X Li C Z Yu D P Liao Z M 2016 Nat. Commun. 7 10769
[42] Yan Y Liao Z M Ke X Van Tendeloo G Wang Q Sun D Yao W Zhou S Zhang L Wu H C Yu D P 2014 Nano Lett. 14 4389
[43] Yan Y Wang L X Ke X Van Tendeloo G Wu X S Yu D P Liao Z M 2014 Sci. Rep. 4 3817
[44] Li C Z Li J G Wang L X Zhang L Zhang J M Yu D Liao Z M 2016 ACS Nano 10 6020
[45] Zhang L Yan Y Wu H C Yu D Liao Z M 2016 ACS Nano 10 3816
[46] Hosur P 2011 Phys. Rev. B 83 035309
[47] McIver J Hsieh D Steinberg H Jarillo-Herrero P Gedik N 2012 Nat. Nanotech. 7 96
[48] Cao Y Waugh J Zhang X Luo J W Wang Q Reber T Mo S Xu Z Yang A Schneeloch J 2013 Nat. Phy. 9 499
[49] Chen C Xie Z Feng Y Yi H Liang A He S Mou D He J Peng Y Liu X 2013 Sci. Rep. 3 2411
[50] Park S R Han J Kim C Koh Y Y Kim C Lee H Choi H J Han J H Lee K D Hur N J 2012 Phys. Rev. Lett. 108 046805
[51] Zhang H Liu C X Zhang S C 2013 Phys. Rev. Lett. 111 066801
[52] Raghu S Chung S B Qi X L Zhang S C 2010 Phys. Rev. Lett. 104 116401
[53] Wang Z Weng H Wu Q Dai X Fang Z 2013 Phys. Rev. B 88 125427
[54] Gorbar E Miransky V Shovkovy I 2013 Phys. Rev. B 88 165105
[55] Kim H J Kim K S Wang J F Sasaki M Satoh N Ohnishi A Kitaura M Yang M Li L 2013 Phys. Rev. Lett. 111 246603
[56] Jeon S Zhou B B Gyenis A Feldman B E Kimchi I Potter A C Gibson Q D Cava R J Vishwanath A Yazdani A 2014 Nat. Mater. 13 851
[57] Xu S-Y Liu C Kushwaha S K Sankar R Krizan J W Belopolski I Neupane M Bian G Alidoust N Chang T R 2015 Science 347 294
[58] Bardarson J H Brouwer P Moore J 2010 Phys. Rev. Lett. 105 156803
[59] Cho S Dellabetta B Zhong R Schneeloch J Liu T Gu G Gilbert M J Mason N 2015 Nat. Commun. 6 7634
[60] Hong S S Zhang Y Cha J J Qi X L Cui Y 2014 Nano Lett. 14 2815
[61] Jauregui L A Pettes M T Rokhinson L P Shi L Chen Y P 2016 Nat. Nanotech. 11 345
[62] Zhang Y Vishwanath A 2010 Phys. Rev. Lett. 105 206601
[63] Washburn S Umbach C Laibowitz R Webb R A 1985 Phys. Rev. B 32 4789
[64] Liu Y Zhang C Yuan X Lei T Wang C Di Sante D Narayan A He L Picozzi S Sanvito S 2015 NPG Asia Mater. 7 e221
[65] Cao J Liang S Zhang C Liu Y Huang J Jin Z Chen Z G Wang Z Wang Q Zhao J 2015 Nat. Commun. 6 7779
[66] Solomon P M Morkoc H 1984 Ieee. T. Electron Dev. 31 1015
[67] Smith R A 1978 Semiconductors Cambridge Cambridge University Press
[68] Du X Tsai S W Maslov D L Hebard A F 2005 Phys. Rev. Lett. 94 166601
[69] Liang T Gibson Q Ali M N Liu M Cava R Ong N 2015 Nat. Mater. 14 280
[70] Narayanan A Watson M Blake S Bruyant N Drigo L Chen Y Prabhakaran D Yan B Felser C Kong T 2015 Phys. Rev. Lett. 114 117201
[71] Feng J Pang Y Wu D Wang Z Weng H Li J Dai X Fang Z Shi Y Lu L 2015 Phys. Rev. B 92 081306
[72] Ishiwata S Shiomi Y Lee J Bahramy M Suzuki T Uchida M Arita R Taguchi Y Tokura Y 2013 Nat. Mater. 12 512
[73] Bangura A Proust C Levallois J Doiron-Leyraud N LeBoeuf D Taillefer L Adachi S Sutherland M Hussey N 2010 Phys. Rev. B 82 020514
[74] Husmann A Betts J Boebinger G Migliori A Rosenbaum T Saboungi M L 2002 Nature 417 421
[75] Geim A K Grigorieva I V 2013 Nature 499 419
[76] Meng J Song H D Li C Z Jin Y Tang L Liu D Liao Z M Xiu F Yu D P 2015 Nanoscale 7 11611
[77] Jin K H Jhi S H 2013 Phys. Rev. B 87 075442
[78] Liu W Peng X Wei X Yang H Stocks G M Zhong J 2013 Phys. Rev. B 87 205315
[79] Zhang J Triola C Rossi E 2014 Phys. Rev. Lett. 112 096802
[80] Amet F Williams J Garcia A Yankowitz M Watanabe K Taniguchi T Goldhaber-Gordon D 2012 Phys. Rev. B 85 073405
[81] Bernevig B A Hughes T L Zhang S C 2006 Science 314 1757
[82] Hsieh D Qian D Wray L Xia Y Hor Y S Cava R J Hasan M Z 2008 Nature 452 970
[83] Kane C L Mele E J 2005 Phys. Rev. Lett. 95 146802
[84] Roth A Brüne C Buhmann H Molenkamp L W Maciejko J Qi X L Zhang S C 2009 Science 325 294
[85] Avsar A Tan J Y Taychatanapat T Balakrishnan J Koon G K Yeo Y Lahiri J Carvalho A Rodin A S O’Farrell E C Eda G Castro Neto A H Ozyilmaz B 2014 Nat. Commun. 5 4875
[86] Balakrishnan J Koon G K W Jaiswal M Neto A C Özyilmaz B 2013 Nat. Phy. 9 284
[87] Qiao Z Yang S A Feng W Tse W K Ding J Yao Y Wang J Niu Q 2010 Phys. Rev. B 82 161414
[88] Qiao Z Ren W Chen H Bellaiche L Zhang Z MacDonald A H Niu Q 2014 Phys. Rev. Lett. 112 116404
[89] Lundeberg M B Folk J A 2014 Science 346 422
[90] Xiao D Liu G B Feng W Xu X Yao W 2012 Phys. Rev. Lett. 108 196802
[91] Mak K F McGill K L Park J M C Euen P L 2014 Science 344 1489
[92] Xiao D Yao W Niu Q 2007 Phys. Rev. Lett. 99 236809
[93] Gorbachev R Song J Yu G Kretinin A Withers F Cao Y Mishchenko A Grigorieva I Novoselov K Levitov L 2014 Science 346 448
[94] Sui M Chen G Ma L Shan W-Y Tian D Watanabe K Taniguchi T Jin X Yao W Xiao D Zhang Y 2015 Nat. Phy. 11 1027
[95] Shimazaki Y Yamamoto M Borzenets I V Watanabe K Taniguchi T Tarucha S 2015 Nat. Phy. 11 1032
[96] Ju L Shi Z Nair N Lv Y Jin C Velasco J Jr. Ojeda-Aristizabal C Bechtel H A Martin M C Zettl A Analytis J Wang F 2015 Nature 520 650
[97] Nielsen H B Ninomiya M 1981 Nucl. Phys. B 185 20
[98] Wan X Turner A M Vishwanath A Savrasov S Y 2011 Phys. Rev. B 83 205101
[99] Weng H M Fang C Fang Z Bernevig B A Dai X 2015 Phy. Rev. X 5 011029
[100] Huang S M Xu S Y Belopolski I Lee C C Chang G Wang B Alidoust N Bian G Neupane M Bansil A 2015 Nat. Commun. 6 7373
[101] Xu S Y Belopolski I Alidoust N Neupane M Bian G Zhang C Sankar R Chang G Yuan Z Lee C C 2015 Science 349 613
[102] Lv B Q Weng H M Fu B B Wang X P Miao H Ma J Richard P Huang X C Zhao L X Chen G F Fang Z Dai X Qian T Ding H 2015 Phys. Rev. X 5 031013
[103] Lv B Q Xu N Weng H M Ma J Z Richard P Huang X C Zhao L X Chen G F Matt C E Bisti F Strocov V N Mesot J Fang Z Dai X Qian T Shi M Ding H 2015 Nat. Phys. 11 724
[104] Xu S Y Alidoust N Belopolski I Yuan Z Bian G Chang T R Zheng H Strocov V N Sanchez D S Chang G Zhang C Mou D Wu Y Huang L Lee C C Huang S M Wang B Bansil A Jeng H T Neupert T Kaminski A Lin H Jia S Hasan M Z 2015 Nat. Phys. 11 748
[105] Yang L X Liu Z K Sun Y Peng H Yang H F Zhang T Zhou B Zhang Y Guo Y F Rahn M Prabhakaran D Hussain Z Mo S K Felser C Yan B Chen Y L 2015 Nat. Phys. 11 728
[106] Huang X Zhao L Long Y Wang P Chen D Yang Z Liang H Xue M Weng H Fang Z Dai X Chen G 2015 Phys. Rev. X 5 031023
[107] Soluyanov A A Gresch D Wang Z Wu Q Troyer M Dai X Bernevig B A 2015 Nature 527 495
[108] Wei P Lee S Lemaitre F Pinel L Cutaia D Cha W Katmis F Zhu Y Heiman D Hone J 2016 Nat. Mater. 15 711
[109] Fu L Kane C L 2008 Phys. Rev. Lett. 100 096407
[110] Das A Ronen Y Most Y Oreg Y Heiblum M Shtrikman H 2012 Nat. Phy. 8 887
[111] Mourik V Zuo K Frolov S M Plissard S Bakkers E Kouwenhoven L 2012 Science 336 1003
[112] Sun H H Zhang K W Hu L H Li C Wang G Y Ma H Y Xu Z A Gao C L Guan D D Li Y Y Liu C Qian D Zhou Y Fu L Li S C Zhang F C Jia J F 2016 Phys. Rev. Lett. 116 257003
[113] Wang H Wang H Liu H Lu H Yang W Jia S Liu X J Xie X Wei J Wang J 2016 Nat. Mater. 15 38
[114] Liu H Neal A T Zhu Z Luo Z Xu X Tománek D Ye P D 2014 ACS Nano 8 4033
[115] Li L Yu Y Ye G J Ge Q Ou X Wu H Feng D Chen X H Zhang Y 2014 Nat. Nanotech. 9 372
[116] Li L Yang F Ye G J Zhang Z Zhu Z Lou W Zhou X Li L Watanabe K Taniguchi T Chang K Wang Y Chen X H Zhang Y 2016 Nat. Nanotech. 11 593
[117] Li L Ye G J Tran V Fei R Chen G Wang H Wang J Watanabe K Taniguchi T Yang L Chen X H Zhang Y 2015 Nat. Nanotech. 10 608
[118] Dean C R Young A F Meric I Lee C Wang L Sorgenfrei S Watanabe K Taniguchi T Kim P Shepard K 2010 Nat. Nanotech. 5 722